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AN EXTENSION OF HEWITT’S INVERSION FORMULA AND ITS 
APPLICATION TO FLUCTUATION THEORY 


§erban E. Badilf0| 


Abstract. We analyze fluctuations of random walks with generally distributed increments. 

Integral representations for key performance measures are obtained by extending an inversion 
theorem of Hewitt hd for Laplace-Stieltjes transforms. Another important step in the 
analysis involves introducing the so-called harmonic measures associated to the walk. It is 
also pointed out that such representations can be explicitly calculated, if one assumes a form 
of rational structure for the increment transform. Applications include, but are not restricted 
to, queueing and insurance risk problems. 
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1. Introduction 

A standard assumption in the theory of fluctuations of random walks as they appear in, 
e.g., queueing or insurance applications, is that the increment of the random walk can be 
represented as the difference of two independent random variables. In the context of queueing 
theory, random walks with increments which do not have this property appear as embedded 
at arrival epochs of customers in a single server queue in which the service requirement of the 
current customer is correlated with the time until the next arrival. In risk reserve processes 
that appear in insurance, the independence assumption is violated when the current claim 
size depends on the time elapsed since the previous arrival, and hence on the premium gained 
meanwhile. 

The purpose of the present paper is to show how much can be done for random walks 
which do not satisfy this independence assumption, regarding their maxima (as in the waiting 
time/maximum aggregate loss), their minima (idle periods/deficit at ruin) and their excur¬ 
sions which are related to busy periods or the time to ruin. 

From a queueing perspective, it turns out that the busy period is a more sensitive issue 
to study than the idle period or the waiting time, as it appears from the proof of Theorem 
[2] below. For this purpose, we present a generalization of an inversion formula for Fourier- 
Stieltjes transforms due to Hewitt HI], which in turn is an extension of P. Levy’s inversion 
formula. All these results are essentially variations on the Dirichlet integral for complex-valued 
functions of bounded variation. 

From a stochastic point of view, the information contained by the increments of the random 
walk is sufficient to infer about the extreme statistics; and since successive increments are 
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independent, the usual form of Hewitt’s inversion formula is sufficient to obtain the integral 
representations; this is how it was used originally by Spitzer m to derive the Laplace-Stieltjes 
transform of the maxima of partial sums. He also related these transforms to the Wiener- 
Hopf problem (see also Cohen [5], Ch. II.5, for the relation with the Wiener-Hopf equation 
as it appears in Probability Theory). For the derivation of the length of an excursion above 
the starting level, there is more information needed, namely, that given by the partial sums 
Y^i=i Bi together with the partial sums of the embedded random walk S n = Bi — Aj. 
It is possible to derive the excursion lengths still using Hewitt’s formula when the A,’s are 
independent of the B^s, and this was carried out in Kingman m • It is shown here that if one 
extends Hewitt’s inversion formula, a similar derivation is possible for the case when there is 
arbitrary correlation inside the vectors (Aj,Bj), which means the random walk S n can have 
generally distributed increments. 

Hewitt’s approach was to find a most general inversion identity for Laplace-Stieltjes trans¬ 
forms with a view on Harmonic Analysis; this is more than needed for our purposes. Instead 
of trying to find a most general instance of inversion, we will focus on obtaining a suffi¬ 
ciently broad result to apply to random walks as they appear for example in the study of 
workload/insurance related problems. One can then hope that the result itself will find ap¬ 
plications in other related areas of probability and statistics. 

In the queueing literature, Conolly [6] obtained (the transform of) the busy period together 
with the number of customers served in a queue with exponential inter-arrivals and indepen¬ 
dent Erlang distributed service times. Conolly’s results from |6j were extended to general 
independent inter-arrivals and service times in Finch [9j and in Kingman m- 

The time to ruin has been studied in the insurance literature by deriving recursion formulae, 
typically obtained by discretising the claim sizes. For example, Dickson and Waters [8] present 
various approximation methods for its numerical computation. Another good reference is 
Prabhu m §3, who obtained an integral equation for the time to ruin starting from a positive 
capital u, in the Cramer-Lundberg risk reserve model; Borovkov and Dickson [4] obtain series 
representations for the distribution of the time to ruin in the Sparre-Andersen risk reserve 
model with exponentially distributed claims and general renewal inter-arrivals. Besides these 
results, there exists a significant amount of literature on the Gerber-Shiu functions which 
contain the time to ruin as a special case. 

The paper is organized in the following way: In Section [2], we extend Hewitt’s inversion 
formula to allow for probability distribution functions on R 2 which do not have a product form 
(see Remark [1]). One of the ingredients of the proof consists of having a precise meaning for 
the conditional distribution of B given A; this is settled as a preliminary. The approach used 
for studying fluctuations of random walks involves obtaining integral representations for the 
above-mentioned quantities. The busy period, idle period, transient workload, and the related 
insurance functionals can still be determined in the form of a Cauchy integral, once Theorem 
|T] is combined with a version of Spitzer-Baxter’s identity (Proposition [[]), and this is carried 
out in detail in Section [3] for the correlated GI/G/1 queue and the Sparre-Andersen risk 
reserve process. Roughly speaking, all the transforms of the relevant performance measures 
are obtained by reading the inversion formula in Theorem |T] from right to left. 

Having obtained a Cauchy integral representation for the busy period, one can then evaluate 
it when the transform of the generic pair ( B,A ) is a rational function in the argument that 
corresponds to the service requirement B (Section dj). 
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2. On Hewitt’s inversion formula 

The starting point is P. Levy’s inversion formula which gives a precise form to the well 
known assertion that a characteristic function uniquely determines a probability measure 4> 
on the real line: 


iT 


lim - 

T->oo 2iri 


-iT 


OO 

J e^ x cj)(da 


e — e 


d £ = ^<K(M)) + 7;<£((<*, &])■ 


Thus the value 4>((a,b]) := f X{a,b\{ u ) du) can be recovered, X(a,b\ being the indicator func¬ 
tion of the interval ( a,b \. Hewitt’s formula extends this result to recover directly functionals 
of the form 


<£(/) : = l J [/(«+) + /(«-)]^(d«), 

for functions f of bounded variation; if f is also continuous, then the integral above becomes 
f f(u)(t>(du). 

The inversion formula in the form given by Levy is further generalized to higher dimensions 
and some other topological groups in Hewitt DU- It is not however the purpose of the current 
paper to explore the possibility of a most general form for it. Such an attempt may not even 
yield satisfactory results, as was already pointed out in m 

We will generalize the above in Theorem |T] to probability measures on M 2 related to the 
random vector ( B,A ). The probabilistic structure is never really lost, the conditional distri¬ 
bution of B given A appears throughout the proof of Theorem |T] disguised as the conditional 
kernel q(u,y), which is defined below as a ’’partial” Radon-Nikodym derivative. 

Preliminaries. We say that / : R —> C is of bounded variation if both its real and imaginary 
parts are of bounded variation; this is the same as |/| being of bounded variation because all 
norms are equivalent on C. We will also work with (complex-valued) measures on the Borel 
subsets of M 2 , which are of finite total variation. For such a measure (j), we will denote by 
\(j)\ the total variation measure of (j>, and with ||</>|| its total variation. This will suffice for our 
purposes, but this set-up is fully detailed and generalized in DU, and the references therein. 

Let (B,A) be a random vector on some probability space, having an arbitrary probability 
distribution. Denote by P the probability measure and by H the joint c.d.f. of ( B,A ): 

H(x,y ) = P (B < x, A < y). 

The correlation device for the increment of the random walk is given by this distribution, not 
necessarily having a product form, and its Fourier-Stieltjes transform 

„(»!„,) := E e-— = »^(dx,d s) . 

In general this is convergent only for IZe s\ = 0, 1Zes2 = 0, but if ( B,A ) is supported on the 
non-negative quadrant in M 2 , then h can be continued analytically to IZe s\ > 0, 1Zes2 > 0. 
The characteristic function of the increment A — B will also be relevant: 




e^dP (A-B<x), 7£e£ = 0. 
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Let A be the probability measure associated with the random vector (B,A), and u be the 
marginal measure associated with B , 


A(WxV) :=P(BGW,deV), u(U) := \{U x R), ZL,VeH(M), 

with B(M) the family of Borel sets on the line. The notation H(du, y ) will be used to suggest 
that we are integrating w.r.t. the measure A y (U) := \(U x (— oo,y]). We will work with a 
version of the conditional cumulative distribution function (c.d.f.) of A given B and this is 
made precise below. 

It clearly holds that A y (U) < v(U), in particular X y <C u, so let 

q(u,y) ■= 

be its Radon-Nikodym derivative. Heuristically, q(u, y) is to be regarded as q(u, y) = P(A < 
y | B £ dri), and we have the disintegration identity 

[ H(du,y) = H(U,y) = [ q(u,y)v(du), 

Ju Ju 

with any of the terms above meaning P(A < y, B £ U). We will be working with a regular 
version of q, which exists by virtue of the separability of R, see for instance Kallenberg m 
Thm. 5.3, p. 84. Furthermore, we have more than just regularity for this kernel, the same 
result gives that q(u,y) is regularly monotone as a function in the argument y, i.e., q(u,y ) is 
non-decreasing in y, for u outside a set of v— measure zero which does not depend on y. 

These considerations are quite intuitive because of the probabilistic nature of the measure 
associated with H. It turns out, however, that we will have to consider instances of the 
inversion theorem for the slightly more general case of complex valued functions H which are 
also of bounded variation, and for this purpose we will show below that H (d«, x) can be given 
a meaning in an analogous way. 

Confusing H with its associated complex-valued measure, see Hewitt HD. we can reduce 
it to a monotonically increasing function, by splitting it into real and imaginary parts and 
using the Jordan decomposition 


TZeH = (UeH) + - (KeH)~, 

for the signed measure TZeH , and similarly for XmH. Setting v±(U) := (JZeH)±(U x M), it 
holds with the similar notation as in the probabilistic case that (7 ZeH)± <C , so we can 
define again the Radon-Nikodym derivatives 


Qi(u,y) ■= 


d (TZeH)+ f ^ ^ d {JZeH)~ 

d„+ W ' 9l (u,y) ~ diF 


(«). 


and similarly for ZmH. Now we can reconstruct H(du,y) in an obvious way, using the 
linearity of the Radon-Nikodym derivative. 

Alternatively, we could have used the total variation measures: 


\1ZeH y \(U) < \KeH\(U x R) 

and thus 7 ZeH y <C u, but a simple argument relying on the Hahn-Jordan decomposition 
shows that this construction yields the same result for H(du,x). 
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Moreover, the monotone regularity property of the probabilistic instance extends to q(u, y ) 
being of bounded variation in y outside a set of u-measure zero which does not depend on y. 
This property will be useful in the proof of the next result. 

Theorem 1 (Hewitt inversion extended). Let H be a totally bounded (complex-valued) 
measure on M 2 , and let f : R —>• C be a function of bounded variation which is also absolutely 
integrable (w.r.t. the Lebesgue measure). Then the following Cauchy principal value can be 
represented as a Lebesgue-Stieltjes integral: 


iT 


Je^ u -^H(du,y)f(y)dy\dC 


-iT 


{—oo — oo 


oo 

\ J {f{u+)H(du,u+)+f(u-)H(du,u-)}. 

— OO 


/ need not be integrable w.r.t. H(du,u). If one of the sides above converges, so does the 
other one. 


Remark 1. If H is of the form H 1 H 2 , the double integral inside the Cauchy principal value 
factorizes into the Fourier-Stieltjes transform of H 1 and the Fourier transform of f(y)H 2 (y), 
this function being again absolutely integrable and of bounded variation. Thus the above 
reduces to the inversion formula in Hewitt ED. Thru. (3.1.1): 


iT f 00 ^ f oo ^ 00 

/ I / e ^ UHl ^ r | / d y I A ^ = \ j ld{u+) + g(u-)\ Hi(du), 

—iT V—00 J V—00 ) —00 

with g(y) = f(y)H 2 (y). 

Proof of Theorem [B For fixed u, change the variable x := u — y, so that the double integral 
inside the Cauchy principal value becomes 


OO OO OOOO 

// e £(w y^H(du,y)f(y)dy = II e^ x H (du, u — x)f(u — x ) dx. 


— OO —OO 

We can bound 


—00 —OO 



H(du, u 


x)f(u 



x)\ |.H|(d-u, 00 ) 


H 


\f(v)\ du, 


hence j H(du,u — x)f(u — x) is absolutely integrable in x because / is. We can now change 
the order of integration in the Cauchy principal value, which becomes after integrating over 

c 


iT 

lim [f — [ e^ x d£,H(du,u — x)f(x)dx= lim [ ^-^- dx [H(du, u — x)f(u — x). (1) 
T— TOO JJ 2iri J T— TOO J nx J 

-iT 

At this point, disintegrate the kernel H(du,u — x) = q(u,u — x)u(du), so that, by the 
Radon-Nikodym Theorem, we can rewrite the right-hand side in dTJ) as 
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lim 

T—>• oo 


[ u(du){ 

[ —- —q{u,u — x)f(u — x) dxl 

1 l 

J ttx J 


(2) 


The main remark is that q(u,u — x)f(u — x ) is of bounded variation as a function in x for 
u-almost all u, and regularity is the key, as can be seen from the following: 

Let {xi}i£i be some ordered sequence determined by the edges of an interval partition of 
R. We can write: 


Y \q(u, u - Xi)f(u - Xi)-q(u, u - x i+ i)f(u - x i+1 )| 
iei 

< Y,\ q ( u ’ u ~ x i)Wf( u ~ x i) ~ f( u ~ x i+ l)l 

iei 

+ Y l/( w - x i+ i)\\q(u,u - Xi) - q(u, u — x i+ i)|, u £ Q, 
i£l 


where, by virtue of the regularity of q, Q is a v —negligible set, outside which q(u, •) is of 
bounded variation. Since Q does not depend on the choice for the sequence {x^}, we can take 
the supremum over all such sequences, and using that /, q are of bounded variation, gives 
that also q(u, u — x)f(u — x) is of bounded variation in x, for u—almost all u. 

We have arrived at the following limit: 


lim 

T—y oo 


/ 


—oo 


sin T x 
nx 


c p(u , x) dx 


;[<£>(?/, 0+) + ¥>(u,0-)], 


for fixed u and ip(u,x) := q(u,u — x)f(u — x). This identity is known as Dirichlet’s integral. 
The integrability condition for the left-hand side that assures the limit exists is that <p(u, •) be 
of bounded variation. As seen from the above, this assumption is only slightly more general 
than Dirichlet’s original condition of monotonicity for <p(u, •). See Doetsch [7| Ch. 24, or 
Titchmarsh [21] 13.2 (the condition is also known as Jordan’s test). 

The proof will be complete as soon as we show that the limit in T can be taken inside the 
u(d u) integral ©. Reduce tp(u, •) to a positive monotonically decreasing function. Using the 
second mean value theorem, for any f3 > 0, we can find j3 > a > 0 such that 


J X 



sin T x 

x 


dx 


7T 

„ x , f sinx , 

< \ip(u, 0) / -dx. 

J x 


0 0 0 
The final upper bound is obtained by a change of variable x —> x/T. The integral along 
the negative half-axis can be treated in a similar way, and this is sufficient to allow the 
interchanging of limit and integration in ©; this together with the fact that Dirichlet’s 
integral identity holds for u-almost all u completes the proof. 


3. The analysis of fluctuations 

Having settled the inversion result in the previous section, let us start with the study of 
the special type of random walk which was briefly described in the introduction. 
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It is assumed that the random variable X that stands for the generic increment of the 
random walk {S n } n >o 


S n = Sq -\- y j Xi , 

i=l 

can be written as the difference B—A, with ( B , A) having some general distribution supported 
on the non-negative quadrant of M 2 , as in Section [2] Moreover, we take the law P of {S^} to 
be the law conditional on So = 0 . Also set b n = Y^i=i an d = Ya= i A'- 

In queueing terms, one can think of the pair ( B , A) as the service time of a generic customer 
together with the time until the arrival of the next customer in a GI/G/1 queueing system - 
which can be always normalized to unit server speed without losing generality. 

Let N be a random variable which is distributed as the number of increments before S n 
becomes negative for the first time. In the language of the single server queue, N is distributed 
as the number of customers served during a busy cycle of the server. Then (assuming unit 
server speed) 6 /v, ajy, —Sn stand for the length of the busy period, that of the busy cycle 
and respectively the idle time of the server. 

In terms of insurance and risk theory, this pair can be interpreted as the time elapsed (and 
hence the premium B gained) since the last claim incurred together with the amount A 
claimed through an insurance policy. Then, conditional on starting with 0 initial capital, 6 jv 
is the time to ruin, a at is the total amount claimed until ruin (including the claim that causes 
it) and —Sn is the deficit at ruin. 

The representation given below was obtained in Wendel [23] §4 as an algebraic identity, 
slightly more general than Spitzer’s identity, who originally obtained in [20] the representation 
for the LST of the successive maxima of the partial sums S n (see also Baxter and Donsker [2] 
for a similar derivation that holds for Levy processes). Theorem 1 and Identity (9) in Kingman 
m are much closer to our purposes. For the sake of completeness, we cite the relevant result 
in the following proposition 

Proposition 1 (Spitzer, Wendel, Baxter). With the above notations, it holds that 

{ °° n I 

- XJ ^- E [ e ~" lbri “ S2S7ll {5™<0}]| » ( 3 ) 

which is valid for Tie si > 0, lie s 2 < 0, \z\ < 1. 

In Kingman [13], this identity is obtained by stopping inside the Spitzer-Wendel identity, i.e. 
in the three-dimensional space where ( a n , b n , S n ) n is evolving, replace the reflecting hyperplane 
at X 3 = 0 with an absorbing hyperplane. This is formally carried out in m by replacing the 
projection operator used by Wendel [23)] with an absorption operator. 

We will use Theorem[j]in conjunction with the version of Spitzer’s identity from Proposition 
[□to obtain integral representations for the transforms of the busy period, idle period and the 
number of customers served during a busy period. 

Throughout the rest of the paper, we will use the dashed integral sign as a replacement for 
the cumbersome limit notation 

iT zoo 

lim / = / . 

T^oo J J 

-iT 


—zoo 
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Theorem 2. We have the following integral representations for P := bjy and I := —Sjv, 
valid whenever Ties > 0, \z\ < 1: 


E{z 


N-sP 


} = 1 — exp 


ZOO 

h / 


s~f 


{log[l - zh(f, 0)] - log[l - zh{f, s - 0]} 


( 4 ) 


E {z N e~ sI } = 1 - exp 


1 

27ri 


j log[l - zh{£, -0] 


( 5 ) 


V —zoo ) 

Here log is the principal branch, that has the cut taken along the negative real axis between 0 
and oo, so that it admits the power series representation 


log 


1 — r 


ry>' v 

M < i. 


n =1 


Proof. Below we will use an integration by parts argument and for this reason it will be 
convenient to introduce the function G(x, y ) = P(-B < x,A > y). We can write by integrating 
over the possible values of b n : 


OO 


OO 


< a , 


H 


G n (dx,x) and 


/ 


< a n ) = / G n (dx,x~), 


o o 

where G n (x,y) := H* n (x, oo) — H* n (x,y), H* n is the n-fold convolution of H with itself and 
G n (dx,y) is the associated integral kernel, as described in Section [2l 


G n (dx,y) = P (b n € dx, a n > y), G n (dx,y~) = P (b n € dx, a n > y). 

Let us begin with (J5J), which means we start with ([3]) for S 2 = 0 and si = s. The first step 
is to represent the expected value inside the series in ([3]): 


-E[e sbn (l{s n <o} + l{s n <o})] — J e sx [-jG n (dx,x —) + -G n (dx,x+)]. (6) 

o 

Before we can use the inversion result from Theorem |T] define the following 2 -harmonic 
measure associated to the sequence G n (dx,dy), n> 1: 

OO 72 

H*(dx,dy) = ^2—G n (dx,dy), (7) 

Z —^ n 

n= 1 

so that, in particular, 



HI is a complex valued measure proper (i.e. it has finite total variation) for \z\ < 1, 



(8) 
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moreover, its LST equals 


/ OO 

e- s ^y H * z (d X , d y) = -J2 h n (si, s 2 ), (9) 

n =1 

since it holds that G n (dx, dy) = —H* n (dx,dy), and where the interchanging of the integral 
with the series is allowed because of absolute integrability: 


0 -s 1 x-s 2 y | 


OO I 177, pp 

H*\(dx,dy) < ^ — // H* n (dx,dy ) = log 

n=l n 



Now let us use Theorem |T] for f(y) = e sy X[o,oo )(y), > 0, which means we can write 

via © and (©: 


Y — E [e _sbn ( 1 {5 ?l <o} + !{ 5 „<o})] 

n=l 


ZOO f OO OO 

^-j\jj^ x ~ [sH)VH *^ x iy) A y 

—zoo v 0 0 > 


d£- 


( 10 ) 


Assume for simplicity that P(A = B) = 0 (i.e. P (S n = 0) is null for all n). This assumption 
is not essential, see for example the discussion in Cohen [5], p. 284. Then the normalized 
indicator function |(l{s„<o} + l{S n <o}) that appears on the left-hand side of ([ 6 ]) simplifies to 
1 {5 7 i <o}. Since H*(dx,y) is of bounded variation in y, we can use the integration by parts 
formula for Lebesgue-Stieltjes integrals, so that (flOl) becomes 


ZOO 

1 f dg 

27ri jT s + £ 

—zoo 



H 


OO 


y =o 


1 f dg 
2m J s + £ 

—zoo 



( 11 ) 


where we used ([9]) and the identity G(x,0—) = F(B < x). Changing the variable £ —> — 
the exponent in (© can be rewritten via (flTT) . for Ties > 0 , \z\ < 1: 


= T j -^{log[l-^(f, 0 )]-log[l-^K.«-£)]}. (12) 

n ~ 1 —io o 

Thus (|4j) follows from these considerations and Spitzer’s identity (JTJ). 

For the integral representation ()5j) . the extension of Hewitt’s formula is not needed. The 
starting point is © with si = 0, S 2 = — s, Tie s > 0, together with the identity 

o+ o- 

Ee s 5 n {5n<0} = l J e M dP (S n <x) + ± j e sx d¥(S n <x ) 

—oo —OO 

valid because P(Si = 0) = 0. Similarly as (J7J), set F{x) = P(5i < x) and introduce 
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00 n 

F:( dx):=J2^F* n (dx). (13) 

n=l 

Use the inversion formula in Remark [T| with g(y) = e sy X(-<x>,o ](y), 7 Zes > 0, so that we can 
write similarly as for (110D - (I11[) : 


200 ( OO 


E f Ee ' S " 1 (Sn<0) Vif\S e& - F;(dl) f / ^ 

n ~ 1 —ioo V—00 ) v 00 

-h 

-L n=l 

200 

= - 77 — / -377 l°g[l — Z H& — £)]> Kes> 0 , 
27 T* J s + 4 


(14) 


after the change of variable £ —> — 

Once ([HI) is replaced into ©, it immediately yields ©• The proof is complete. 

Remark 2. If B is independent of A. then h(s 1 , S 2 ) is the product of the marginal transforms 
and the integral representation ([1]) reduces to that from Kingman P3!> Thm. 4 (see also 
Cohen [5], p. 304 for ©). 

Remark 3. The integral representations (14|) . d5]) hold under very general conditions (there are 
no regularity assumptions for the distribution of (B,A), these can even be discrete random 
variables, in which case the LSTs become generating functions). The reason is that these 
representations are given for the interior of their convergence domains (7 Zes > 0, \z\ < 1). 
If we want to take any of the arguments to their respective boundary, we have to require 
extra conditions to ensure convergence. For example, when letting s converge towards the 
imaginary axis, there is a singularity appearing, because the factor l/(£ — s ) gains a simple 
pole located at £ = s. 

It turns out one can give a definite meaning to these integrals, for 7 Ze s = 0, if they are 
regarded as singular integrals w.r.t. the Cauchy kernel l/(£ — s). Then one considers the 
Cauchy principal value obtained by removing a circle of arbitrarily small radius around the 
singularity s and then taking its radius to 0. Now we are dealing with a double principal value: 
first one coming from the pole at s of the Riemann integral along the segment L := [— iT,iT] 
(T large enough, so that s £ L) and the second one obtained by letting T —>• 00 . A standard 
condition (see Gakhov [TO], Muskhelishvili m) that ensures the first principal value converges 
is that the density functions Vh(0 := log[l — 7i(£, — £)], ^(s) := log[l — h(£, s — £)] are Holder 
continuous along the imaginary axis, with some positive indices. In the case of (|26l) . the Holder 
continuity of </q(£) is fairly close to Spitzer’s [19] integrability condition, which requires (upon 
taking s —> 0 ) that [1 —/i(£, —£)]/£ be integrable on a neighbourhood of s = 0 on the imaginary 


axis. 
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The number of arrivals during an excursion. Further along the lines of Remark [3l we 
will use the doubly dashed integral sign to denote the double Cauchy principal value. The 
choice for the branch of log is essential for the definiteness of the first principal value. It turns 
out to be convenient to work with a branch that has the cut between 0 and oo taken inside 
the negative half-plane, so we have by definition, 


_L f ^... _!*.> + dd |log(!r _- logHr -. )] + ±/ (18) 

L L 

The integral on the right is well defined as a Riemann integral as soon as <p(s) is Holder 
continuous along the line L. By choice of logarithm, the argument of log(fT — s) — log(— iT — s) 
equals iri for all T. This means that the first two terms cancel in the limit T —>• oo, and the 
above definition becomes 


ZOO 



—zoo 


zoo 

1 f <p(Q - 
2iri J £ — s 

—zoo 


(16) 


Formula (|15jl differs from the definition given in Mushkelishvili |T5], P- 27 or in Gakhov 
[TO] . p. 16, because therein the cut of the logarithm is taken in the opposite half-plane. To 
be more precise, in our case, Cauchy’s integral representation reads: 


ZOO 



—zoo 


2iri, s £ L + , 
0, s £ L , 

0, s £ L. 


The first two values are the well known Cauchy integral identities; the third one is the 
Cauchy principal value for this specific choice of log (using (fl6l) with ip = 1). 

It helps to think about the Riemann sphere as the one point compactification of the complex 
plane so that the imaginary axis is closed into a large circle on the sphere. Then all of the 
conventions above are specifying a means of integrating on the large circle of the sphere, the 
integrands being extended by continuity at infinity; the interior of the imaginary axis is by 
definition the left hemisphere and the exterior is the right hemisphere. Define 


= j a € C, 

—zoo 

and the integral is defined in the sense of (fT6l) . when Ties = 0. Let us denote by d> _ (s) 
the limit as s approaches the imaginary axis from its exterior (the right half-plane), and by 
$+(a), the limit taken from the interior. Then the Plemelj-Sokhotski formulae (cf. Gakhov 
[10] . p. 25) become with the above conventions: 

<h~(s) = <3?(s), <h + (,s) = <3?(s) + ip(s), 7les = 0. (17) 

Now we can calculate the limit as s approaches the imaginary axis from its exterior in (|4]), and 
from its interior in (J5]h using (1171) . Still denoting the interior limit with l> + (s), © becomes 
for |^| < 1 and as s tends to the imaginary axis: 
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<h + (s) = 1 — [1 — zh(—s, s)] exp 



d£ 

e + s 


log[l - zh(£, -0] 


(18) 


V —200 / 

Before we give the limiting values of ([4]) . we point out how one can simplify it. Consider the 
first integral in the exponent of dH): 


^C 7 ?) := _ 27ri / _ ^(^°)]’ 7Zer l >0 - ( 19 ) 

—200 

The key remark is that h(£, 0) = Ee - ^ is analytic for IZef > 0 and this implies that the 
density log[l — zh(£, 0)] is again analytic for 7 Ze£ > 0, because of the choice of log. Since the 
integrand has a simple pole in the positive half-plane located at rj = £, it follows at once from 
Cauchy’s Theorem applied to the imaginary axis that T (r/) equals 


T(r?) = log[l - zh{r],0)}, 

for any IZer) > 0. This can be seen by closing the segment [— iT,iT ] with the half circle 
spanning between its endpoints inside the positive half-plane. The contour integral thus 
obtained equals log[l — 2/i(s,0)] for T large enough, so that the pole rj = £ lies inside the 
contour (remark that because of the conventions on the interior of the imaginary axis, this 
contour is traversed in the clockwise direction). Finally, the contribution along the half-circle 
tends to 0 as T —> oo because the integrand behaves as o(|^| _1 ) along the half-circle. 

Having settled the first term in the exponent of Q, we can use the Plemelj-Sokhotski’s 
formula (the continuity of the exterior limit from (|17D ) for the other term, to obtain the 
identity 


$ (s) = 1 - exp | log[l - zh(s, 0)] + 

V —200 

hence we can rewrite (1201) as 


Z-s 


log[l - zh(£,s-£)}}, lies = 0, (20) 


{ 200 1 

7 ^- j 1 !og[l - zh(X, s - 0] > , Hes = 0. (21) 

—200 / 

In particular, for s = 0, the limits (I18j) and (121 j) agree and these must then coincide with 
the generating function of N. The following has been proven 


Proposition 2. With the above notations and conventions, it holds for \z\ < 1: 

d£ 1 -zh(£,-£) 


E z N = 1 — (1 — 


n 'i J 1 f d £ l-zi 

(1_ " )exp ss / t g —r 


— z 


( 22 ) 
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One can determine from Spitzer-Baxter’s identity (see Spitzer m, Thm. 3.1) the transform 
of the transient waiting time in a similar way as above. The Spitzer-Baxter identity reads for 
M n = maxl^o, S u ..., S n }: 

{ o° n ^ 

—^e~ sSt \ , (23) 

n= 1 J 

with 5+ = max(S , n , 0). Keeping the assumption that P(5'i = 0) = 0, we can write 

oo « oo 

Ee- sS " = / e- sx X[ 0,o O )(x)dF;(x) + Yl 

n=l ' n =1 

for F* defined in (fT3l) . It is worth to point out that in this case 5+ has the mass of S n 
distributed along (—oo,0] swept into the origin as an atom. The indicator function appearing 
in the last term in (124[i is not absolutely integrable w.r.t. the Lebesgue measure, therefore 
introduce g e {x) = e £X X(-oo,o){ x ), e > 0, in order to be able to use Hewitt’s inversion formula. 
This same perturbation of the indicator function was used by Spitzer [20]. We have similarly 
as (HU): 



z n C 

— J X(-oc,o){x)dF(S n <x). (24) 



OO yi -i /» l v- 

n ^ —io o 

We can take e —> 0, by using the Plenrelj-Sokhotski identity (1171) for the right-hand side and 
dominated convergence on the left ( \z\ < 1): 


v — ns n < o) 

n 


n= 1 


2m J £ 1 — z 


(25) 


Hewitt’s inversion formula is directly applicable for the first term in (1241) : as in the proof of 
d5]) , use Remark [T| for g(y ) = e~ sy X[o,co)(y)i Fes > 0: 



200 



— ZOO 


__ OO 

OO y) p 

Ee^» 

n= 1 11 n 


e -(o+Cydy i d£ 


ZOO 



—200 



after the change of variable £ This together with (1231) yields 


j^jz n Ee 


—sM n 


n =0 


1-Z 


exp 




V —zoo —zoo 


log 


1 ~zh(£,-£ ) 
1-z 


(26) 

with Fes > 0, \z\ < 1 (cf. Cohen [5], (5.29) p. 276). This formula is related to that obtained 
by Spitzer |2D] concerning the Wiener-Hopf equation which has a probability density as the 
kernel. 
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4. Examples 

In this section we evaluate the integral representation (j4|) . under the assumption that the 
transform of the generic pair ( B , A) is a rational function in the argument that corresponds 
to the service requirement B. The analysis is similar to the one carried out in Cohen [5], Ch. 

II.5. 

Assume that for all s 2 , the joint LST /i(si,s 2 ) is a rational function in the argument si, 
which can be represented as 


h(si,s 2 ) 


hi(si,s 2 ) 

h2(si,S2) 1 


(27) 


where /i*(-, s 2 ) are polynomial functions. Moreover we will assume that for IZe s >o, M&s-O 
has a finite number of poles in the negative half-plane as a function in the argument £ 
(/i(£, s — £) is already meromorphic in this region, because of the previous assumption). This 
is not essential, but an algorithmically friendly assumption which will give a representation 
for the busy period transform in terms of a finite number of factors. 

We may still assume without losing generality that P(.B — A = 0) = 0, which implies 
h(s\,S 2 ) —>• 0, as si —>• oo, !Zes\ > 0, and the convergence is uniform in s 2 , for 1Zes2 > 0. In 
particular, we have for any TZes 2 > 0, deg/ii(-, s 2 ) < deg /i 2 (-, s 2 ). 

Before we proceed with the analysis, let us point out some ways of creating correlation 
between the inter-arrivals and the corresponding service times. 


Example 1 (Threshold dependence) This is one of the simplest ways of making B depend 
on the size of A: for a fixed threshold l > 0, B ~ B\ on the event A < l and B ~ B 2 
otherwise; with Bi independent of A and having rational transforms fi(s 1 ), * = 1,2; thus 


h{s 1 ,s 2 ) = /i(si)ai(s 2 ) + f 2 (si)a 2 {s 2 ), 

l OO 

ai(s 2 ) = f e~ S2V dF(A<y), a 2 (s 2 ) = f e~ S2V dF(A<y), 
y =0 y=l+ 

so that ai(si) is an entire function and a 2 (s 2 ) is analytic and bounded for lZes 2 > 0. This 
construction can be naturally extended to k thresholds, giving 

k 

h(si,s 2 ) = ^2 fi(si)ai(s 2 ), 

Z— 1 

with aj(s 2 ) entire functions, i < k, and ak{s 2 ) is analytic for TZe s 2 > 0. 

Example 2 (Markov Modulation) Let (X n ) be a finite state Markov chain which has an 
absorbing state and denote with k the number of jumps until absorbtion. Define 


K, 

(A,B) = J2(Ai,Bi), 

i= 1 

where (A,;,Bj) are i.i.d. vectors. The component A± is allowed to be generally distributed 
with 5 o( s 2 ) = Ee -S2j41 , and B\ has a rational transform of the form /i(si)// 2 (si). 
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If we denote by ot the initial distribution of (X n ), by T the transient component of its 
transition matrix, and by t the vector of exit probabilities, then by conditioning on At, the 
transform of ( A,B ) is (see for instance [3]): 


h(si,s 2 ) = of 


./Ml) 


i -i 


Ji{si)g 0 (s 2 ) 

Both these examples are of the form assumed by (1271) . 


I-T 


t. 


Remark that h 2 (-,s 2 ) can only have zeroes with negative real part, because of the regularity 
domain of h(-,s 2 ). With these assumptions, the exponent in Q becomes 


ZOO zoo 

log[l-IE^e-* p ] = ^- f -^log[l-^(£, 0 )]-^ -f ^-log[l -zh^s-O}- (28) 

—zoo —zoo 

For the principal branch of the logarithm which has the cut taken along the negative real 
axis, the single valued functions log[l — z/r(£, 0 )] and log[l — zh(£,s — £)] are holomorphic, 
for £ lying in a neighbourhood of infinity, 1Ze£ < 0. The reason is that for such values of 
£, |z/i(£, 0)| < 1 , | zh(£,s — £)| < 1 , and then a simple geometric argument shows that both 
1 — zh(£, 0) and 1 — zh (£, s — £) lie in the positive half-plane. With this choice of the cut, the 
evaluation of the integrals (1281) becomes an application of the theorem of residues. Before we 
can evaluate ()28[) . the zeroes and poles of the arguments of the logarithm must be localized. 
The following lemma will also be useful later on. 

Lemma 1. The functions h 2 (£, s — £) and h 2 (£, s — £) — zh\(£, s — £) have the same number 
n = n(s,z) of zeroes in the negative half of the complex £,-plane, either when \z\ < 1, IZes > 0, 
or \z\ < 1, Ties > 0. 

Assuming Ties > 0 and KB, EA < oo, then under the extra ergodicity condition E B < E A, 
the functions h 2 (£, s — £) and h 2 (f, s — £) — h\(£, s — £) have the same number m = m(s) of 
zeroes with negative real part. 

Proof. The statements will follow from Rouche’s theorem (cf. Titchmarsh [21] . p. 116) as 
soon as we show that these functions are analytic in the interior of some suitably chosen 
contours and on their boundary it holds that 

M£, s - 01 > a - 01, MO -01 > IM0 -01- (29) 

Fix R > 0 and consider the contour C consisting of the segment of the imaginary axis 
between —iR and iR together with the semicircle with radius R that spans in the negative 
half-plane. For the segment of the imaginary axis, we have the following bounds 

\zh(£,s- 01 = M|Ee- ?p - (s -M < |z|eM b - (s " ?)a | < I^Ee^M Tle^ = 0. (30) 

For the bound on the half-circle, consider the following representation for h(s\,s 2 ): 

w v _ CL\(s 2 )sf 1 + a 2 (s 2 )s'i 2 + ... + a n (s 2 ) 
bi(s 2 )si + b 2 (s 2 )si 1 ... + b n -\-i(s 2 ) 

where n = n(s 2 ) (remark that for the examples presented above, the denominator does 
not depend on s 2 , hence neither does the degree n). The functions ai(s 2 ) can be taken to 
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be bounded for lZes 2 > 0, because it holds that f/i(l,S 2 )| o 0 as S 2 o oo. Moreover, 
we can assume 61 (^ 2 ) = 1, for lZes 2 > 0, after normalizing the fraction. For fixed S 2 , 
let £i(s 2 ) be the zeroes (all with negative real part) of ^ 2 (^ 1 , S 2 ); when bounding the above 
representation of h(s 1 , s 2 ), use the triangle inequality for the numerator and use the inequality 
|— ^ 21 > | |^i | — |^ 2 || for the denominator: 


IMS1,S2)| = 


&(-»)) 


i= 1 


> 


ifm-i&md 


2=1 


the right-hand side being a polynomial function in |si| of the same degree as h, 2 (s 1 , S 2 )- Thus 
we have the upper bound 


Ifr(£,s-OI < 


ai( S -oiier 1 + ia2( S -oiier 2 + ... 


im=i(iei-i&(--oi)i 

Then it follows from the facts that deg/ii(£, s — £) < deg/i 2 (£,s — £) and that the a, are 
bounded, that 


\h(C,s-0\=om- 1 ), \H\ = R^oo,Kes>0, (31) 

for £ running along the half-circle that closes the contour C. 

From the bounds (1301) and (13T1) . it follows that s — £)| < |/i 2 (£, s — £)| when £ is on 

C, for R large enough, either if \z\ < 1, 7 Zes > 0, or if \z\ < 1, IZes > 0. This yields the first 
part of the lemma, via Rouche’s theorem. 

For the second part, consider the contour C e made up from the segment that runs in parallel 
to the imaginary axis and lying to its left at distance e, together with the arc of the circle 
with radius R spanning in the negative half-plane between the edges of this segment. 

It is essential that /&(£,—£) is meromorphic in 7 Ze£ < 0 (see the discussion below ([271) 1. 
Since it has isolated poles, we can find e > 0, such that /i(£, —£) is holomorphic in the thin 
strip —2e < Tie £ < 0. Then the left derivative of the function 7i(£, —£) exists at 0 and we 
have by hypothesis, 

lim 4rh(£, -0 = EA-EB > 0, 
d£ 

Tlet<0 

in particular, h(He£, —7£e£) < h( 0,0) = 1. We can now bound for Res > 0 and £ lying on 
the segment of C e : 


\h(C,s~0\ < E|e" ?B - M+5A | < h(Ke£,-1let;) < 1. 

The bound for £ lying on the arc component of C e follows in the same way as (EB. By virtue 
of Rouche’s theorem, the proof is complete. 

Remark 4. It follows in a similar way as for the first part of Lemma Q] that the polynomials 
7 i 2 (£; 0 ) and /i 2 (£, 0 ) — zhi(£, 0 ) have the same number of zeroes with negative real part, 
\z\ < 1. But since 7 i 2 (£> 0 ) has only such zeroes and deg/ 12 (•,()) > deg7ii(-,0), the same holds 
for h 2 (£, 0 ) - zhi(£,0). 

The idea for evaluating (1281) is to use the theorem of residues for the contour integrals along 
C e while arguing that the contributions from the integrals along the half-circle vanish as the 
radius R —> 00 . Focus on the contour integral of the second term in ([28|) taken along the 
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semi-circle component, say S e , of C e . For R large enough S e will be contained in the interior 
of a domain where h(C,s — £) is holomorphic, and in addition, \zh(C,s — £)| < 1, hence the 
position vector 1 — zh{C,s — £) has positive real part when the argument £ runs along S e . 
This means log[l — zh(C,s — £)] is holomorphic in a neighbourhood around the arc S e . In 
conclusion, we can integrate by parts: 


L 


dtj 


log[1 - zh(£, s - f)] = - log(s - 0 log[l - zh(C, s - 0] 


— e+iR 


—e—iR 


+ I log(s - £)-T7 - zh (C, a - 0]d£- 

JS e 


For large R, \h(C, s — £)| —> 0, which means | log[l — zh(C, s — £)]| ~ | zh(C, s — £)|, so the first 
term on the right behaves in absolute value as 


~ |log(£ - s )| | zh{£,s- 01 ~ \z log(C - s)l 10 \ 
and the integrand on the left-hand side behaves as |£|~ 2 . Thus we have 

/ log(£ - s)x 7 lo g[l - zh(C,s - £)]<!£->• 0, |£| -> oo, UeC, < 0. (32) 

JS e ^ 

Now we are ready to calculate the contour integrals (1281) . Remark that the first term is 
the same as (fT9l) and thus equals log[l — z/i(s,0)], as was found in Section [3l For the second 
integral in (1251) , fix IZe s > 0 and consider the integral taken along the contour C e described 
in the proof of Lemma [TJ e is taken sufficiently small such that no poles of the integrands are 
lying between the segment and the imaginary axis, irrespective of R (this can be found since 
there are finitely many poles in the negative half-plane). The integral can be approximated 
from the interior of the negative half-plane using the contours C e , for arbitrarily large R 
and small e. Splitting the integral along C e based on the factors inside the logarithm and 
integrating by parts in each term, the expression in (I28j) becomes 


1 

2iri 



log(s - C) 


^[h 2 (C,0)-zh 1 {C,0)} 

M6°) - zM£,o) 




ih 2 (C,0) 

0 ) 


d£ 




s - 0 


dC 



h 2 (C,s- 0 - zhiiC, s — C) 


By (1321) . the total contribution from the integral along the semi-circle S e vanishes as R —>• oo. 
Moreover, the branch of log was chosen such that the factors log(s — £) are analytic for 
lZe£, < 0. Then the integrands have simple poles located at the zeroes of their denominators 
in the negative half of the complex plane. Denote by £j(s), Ci(z,s), i = 1, the zeroes 
with negative real part of h 2 (C, s — £), respectively h 2 (C, s — Q — zh\{C, s — £), as functions of 
the variable £ (see Lemma [1]), and with r)j, r]j(z), j = 1, the zeroes (all having negative 

real part) of h 2 (C,0), respectively h 2 (C,0) — zhi(C,0). Then the integral in (1251) equals 


n;=i[s - vj(z)} n: ,> - cm _ [ h2 ( S , 0 ) - ^m)] mut* - cm 

i\7=iis - vj] nr=i- &(*, s)] og h 2 ( s , 0 ) no - &(*, *)] 
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after letting R —> oo and e —> 0. Remark also that the degree n is a (piecewise constant) 
function of the argument s (Lemma [I]). In conclusion, we have 


¥j{z n e sP } = 1 — [1 — z/i(s,0)] 


nti[»-&w] 


'ntii*-&(*.«)]' <33) 

The calculations that led to (1551) can be repeated for z = 1 and IZe s > 0. The contour of 
integration is the same as C e , and the second part of Lemma Q] must be used to conclude about 
the number of zeroes under the condition E A > ELL This condition is necessary for stability 
(P has a proper probability distribution, which must be verified when taking z = 1, s = 0). 
The conclusion is that we are allowed to formally replace z = 1 in ()33[) . which becomes 


Ee _sP = 1 — [1 — /i(s,0)] 


rnijs-&(*)] 


m iis-ah s)v 

with the remark that P has indeed a proper probability distribution. 

Similarly, the exponent in (1221) . which also appears in ()26[) . can be rewritten as 


2m J £ 


log 


1 - z/i(£,-£) 


1-z 

whereas the other integral that appears in 


2 s/ log( ‘ ?) 


becomes 


Hfl 1 -zh(£,-S)] 


zoo 

2 h / 

—zoo 


1 

27 xi 



1 -zh(f;,-£) 


Using the assumption (1271) . it is easy to see that the denominator of both integrands above is of 
the form / 12 (C) — C)[^-2(C) — C) — 2/11 (£> — £)]• A similar analysis as for (1551) yields the generating 
function of the sequence {Ee _sMri } n and the generating function of N from Proposition [2} 


(l-z)J]z n Ee 


—sM n 


E / = 1 - (1 - z) 


n? =1 («-fc(o)) n-=iC,(^o) 
"n^-^, 0 )) n;=iO(o) 1 

M < 1, Kes> 0. 


(34) 


Cj(s), £j(z, s ) are defined as in (1551) and n = n(s, z) is the number of zeroes with negative real 
part, as given by Lemma [T] Under the stability condition E B < E A, lim,,.-^ M n = M exists 
in distribution and we can take z —> 1 in (134[) by virtue of Abel’s theorem (Titchmarsh |21j . 
1 . 22 ): 


w P ~ sM = Cj(o)) rij=iCi(i) 0 ) 

' n;u(s-o(i,o)) n;=iO(o)' 

Finally, it is possible to derive the transform of the idle period using a similar contour 
integral; and the methods of this section equally apply to the symmetric case when the 
transform h(s\,S 2 ) is a rational function in the argument s 2 , for each fixed si, IZesi > 0. 
The analysis relies in this case on localizing the poles and zeroes lying inside the positive-half 
plane. 

Concluding remarks The methods used in Sections [5] and [5] have been first applied 
in Queueing theory by Pollaczeck [10J and many others have followed (see also the survey 
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of Takacs [22]). With respect to calculating the quantities related to the fluctuations of 
random walks, from a general perspective, the analysis presented in this paper shows the 
effectiveness of the Radon measure and its associated integral (what we called the Lebesgue- 
Stieltjes integral) when combined with Spitzer’s identity d3|), via the harmonic measures (J7|) 
and (fl3l) . 

The following functions: 

ipf (s) := (1 — z) ^2 z n Ke~ sMn , 7 Ze s > 0, \z\ < 1, 

and 

( s ) := 1 — Kz N e sI , Ties < 0, \z\ < 1, 

make up the solution of the homogeneous Hilbert problem associated with the Wiener-Hopf 
equation. To be more precise, i(>£(s) and ^j(s) satisfy the boundary relation: 

'ip-(s) ipt (s) = 1 — zh(s, — s), 7 Ze s = 0, \z\ < 1, 

and are unique with the additional property that tpf (s ) is analytic in 7 Ze s > 0 and i/>~ (s) 
is analytic in Ties < 0, both continuous up to the boundary (the imaginary axis). The 
uniqueness holds because the kernel 1 — zh(s, — s) has index zero along the imaginary axis 
(i.e., its logarithm is single-valued along the imaginary axis, see Gakhov [ID] . Ch. II), a fact 
which was used throughout Sections [3] and |H It is easy to check that the boundary relation 
holds, if we use the current version of the Plemelj-Sokhotski identities (1171) . and take s to the 
imaginary axis both in ([5]) and (1261) . 

There is a large amount of literature on Wiener-Hopf equations both in Analysis and 
Probability. A reference in analysis is Krein’s manuscript p3j; for probabilistic background 
and applications, see Asmussen p]. It seems Rapoport [Mj was the first to observe the 
connection between the Wiener-Hopf equation and the theory of Riemann/Hilbert boundary- 
value problems. 
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